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The perturbation theory expansion presented earlier to describe the phase-ordering kinetics in the 
case of a nonconserved scalar order parameter is generalized to the case of the n-vector model. At 
lowest order in this expansion, as in the scalar case, one obtains the theory due to Ohta, Jasnow and 
Kawasaki (OJK). The second-order corrections for the nonequilibrium exponents are worked out 
explicitly in d dimensions and as a function of the number of components n of the order parameter. 
In the formulation developed here the corrections to the OJK results are found to go to zero in the 
large n and d limits. Indeed, the large-d convergence is exponential. 



I. INTRODUCTION 



A perturbation theory expansion for treating the scaling features of phase-ordering kinetics in unstable systemsEl 
was presented earlierB for the case of a nonconserved scalar order parameter. In this paper this method is extended 
to the case of a vector order parameter with n components. It is found that spin-wave degrees of freedom introduce 
some new elements into the theory. However, as in the scalar case, the Ohta, Jasnow and Kawasaki (OJK) theorjoEl 
emerges as the zeroth-order approximation and at this order the nonequilibrium indices A and v do not depend on the 
number of components, n, of the order parameter. The second-order corrections to these exponents are determined 

?j ' as functions of d and n, and these corrections vanish for both large d and n. _ 

The theory developed in I and here is a two-step process which builds on earlier worka. First one maps the original 
problem for the order parameter -0 onto one for an auxiliary field m. The field configurations associated with the 
properly chosen auxiliary field are smoother and therefore easier to treat than for the order parameter. The second 
•<H- step is to treat the nonlinear field theory satisfied by fa. It is found in both scalar and continuous (n > 1) cases 

that one must construct, as one constructs a fixed-point hamiltonian in critical phenomena, the equation of motion 
I" satisfied by the auxiliary field in the scaling regime. A new aspect in the development is the conjecture that there is a 

general relation v = 2A — d connecting the nonequilibrium indices. This relation, and the self-consistent maintenance 
of the t 1 ' 2 growth law for this problem, fixes the form of the auxiliary field equation of motion at second order in the 
perturbation theory. An important consequence of this procedure is that the second-order contributions to the indices 
are exponentially small in the large-<i limit. In the continuous case it is only the longitudinal part of the auxiliary field 

C^ ' equation of motion which must be determined self-consistently. The transverse part is unambiguously determined by 

C^ | consideration of the spin-wave degrees of freedom. 

I ■ The nature of the perturbation expansion introduced in I is elucidated further here. For general n it can be seen 

that one can develop an expansion of the nonlinear terms in the auxiliary field equations of motion in terms of a set of 
vertices labelled by its number of spin labelsEl. A vertex with £ labels can be self-consistently taken to be of0(£/2 — 1). 
This, in turn, leads to the result that the I th order cumulant is also of order 0(£/2 — 1). We will refer to this expansion 
as the vertex expansion. As in the scalar case, this expansion is well behaved in the lowest orders of perturbation 
theory. We obtain nontrivial values for the nonequilibrium exponents at second order in perturbation theory by 
exponentiating diverging logarithms which are driven, for general dimensionality d, by internal time integrations. 



II. OVERVIEW 



We study here the phase-ordering kinetics generated by the time-dependent Ginzburg-Landau (TDGL) model 
itisfied by a nonconserved vector order parameter ibjfr.t), (i = 1, 2 . . . ,ri), 



satisfied by a nonconserved vector order parameter i\)i (r , t) , (i 

d^ = _ T SF 
dt 5ip 






where T is a kinetic coefficient, F is a Ginzburg-Landau effective free energy assumed to be of the form 

n d 



\ i=l or=l 



(2) 



where c > and the potential V is assumed to be of the symmetric degenerate wine-bottle form, V = V(4> 2 )- We 
expect only these general properties of V will be important in our analysis, ff is a thermal noise which is related to 
T by a fluctuation-dissipation theorem. We assume that the quench is from a high temperature (Tj 5h2c), where the 
system is disordered, to zero temperature where the noise can be set to zero (ff — 0). It is believedtl that our final 
results are independent of the exact nature of the initial state, provided it is a disordered state with short-ranged 
correlations. 

If we rescale length and times we can put our equation of motion in the dimensionless form 

A(l)VMl) = ~VS(1)} = - ^ff , (3) 

where the diffusion operator 

A(D = £-V? (4) 

is introduced along with the short-hand notation that 1 denotes (ri,ii). The standard ip 4 form for the potential is 
given by V = -^tf 2 + \(if 2 ) 2 ■ 

For late times, following a quench from the disordered to the ordered phase, the order-parameter dynamics obey 
scalingH governed by a single growing length, L(t), which is characteristic of the spacing between defects. In this 
scaling regime the order-parameter correlation function has a universal scaling form 

C(12) = ($1) • $2)) = ^(x,h/t 2 ) 

where ipo is the magnitude of the order-parameter in the ordered phase. The scaled lengthy is defined as x — 
(ri—r2)/L(T) where, for the non-conserved order parameter case considered here, the gepwth lawO goes as L(T) ~ T 1 / 2 
where T = \(t\ + £2)- ln the case of the autocorrelation function fi = ?% = r we haveQ 



(^txj.^f.ta^wf-^l (5) 

where A is a nontrivial nonequilibrium exponent which enters when either t\ or t 2 is much larger than the other. At 
equal-times the scaling function F(x) = J-(x, I) is a nonanalytic function of x for small |a;|. This is best reflected 
in the generalized form of Porod's lawcllj as expressed in terms of the Fourier transform T(Q) ~ Q~( n+d > for large 
scaled wavenumber Q. The large- a; behavior of the scaling function can, with proper definition of x, be put in the 
form T(x) ~ -^e~^ x where v is a nontrivial subdominant exponent introduced in Ref.( ^). 

Values for A and v were found in I which suggested that they are independent. In earlier work in Refs.( |Ic||ll] ), 
within the approximation developed in Refs.( ^), the two exponents were related by 

v = 2\-d . (6) 

We conjecture here that the relation given by Eq.dq) holds more generally and the perturbation theory must be 
constructed to respect this result. The argument leading to this conjecture follows from the assumption that in the 
bulk scaling regime and away from defects the only available vector is ip(l) and the potential contribution to the order 
parameter equation of motion must be of the form 

vS(l))=- J ^- ) i>i(l) (7) 

where a is an unknown constant and, if we are to have scaling, the coefficient in Eq.ufl) must be proportional to L~ 2 . 
In Refs.( ||) and ( [nj) it was found that a = n/2. With these assumptions one can follow the development in Ref.( [|) 
to show that the exponent v can be written in the form 

v = d-- , (8) 

A* 



where /i is the eigenvalue which plays a prominent part in the theory developed in Ref.( ||). Similarly, generalizing 
the analysis in Ref.( |l0| ) for two-time correlation functions in the limit of large time-separations, one can show that 
the nonequilibrium exponent A is given by 

X-*-± . (9) 

Eliminating j- between Eqs.(||) and (||) we obtain the scaling relation given by Eq.(g). 

III. AUXILIARY FIELD METHOD FOR THE N- VECTOR MODEL 

Our goal here is to generalize the perturbation theory expansion approach developed in I for a scalar order parameter 
to the case of general n. We will again express the order parameter as a sum of two pieces 

if = a + u (10) 

where, as in Ref. O, a is the ordering component of the order parameter which depends locally on an ri-component 
field to: a = a(rh). The field u represents the fluctuations about the ordered configuration. The functional dependence 
of <? on to is determined as a solution of the Euler-Lagrange equation for the associated stationary defect problem 

£^ = Wm)] , (ID 

where to is taken to be the coordinate and the solution must satisfy the boundary condition limi^^^a 2 — ipQ. If we 
introduce the notation 

a i;jiJ2---je 



drrij 1 drrij 2 . . . dmj i 



then Eq.(|ll|) reads 

5> 4; « = ^Hto]] • (12) 

f. 

One can obtain the defect profile analytically for the case of a general degenerate potential for the scalar order 
parameter case. In the particular case of a ip 4 potential, one obtains the usual interfacial kink solution a[m] = 
tanh(m / 'V2) ■ For systems with continuous symmetry, n > 1, one does not have a closed form solution for the defect 
profile even for the ?/> 4 potential. However one can make some general statements about the form of the profile 
in the ordered bulk regime. For the lowest charge defects we can write the order-parameter profile in the form 
a[m\ — A(m)m, then the Euler-Lagrange equation reduces to an equation for the amplitude A(m) given by 

„// (n-1) / „ A\ dV(A) . . 

A" + ^ 1 (A'--)= — y- . (13) 

TO \ TO/ OA 

Solutions of Eq.(|l^) for large to, in the bulk away from defect cores, in contrast to the scalar case which has an 
exponential approach to the ordered value, show an algebraic approach to the ordered value: 

A = ^ (l - j- + . . .J (14) 

where £ 2 = (n — 1)/V"(ijjo)- This solution requires that the ordered value of the magnitude of the order parameter 
be given by the solution to V'(i/jo) = 0; an d, for the solution to be stable, we require V"(ij)o) = q^ > 0. There are 
some additional general properties we need in the large to regime. If we write 

and use the expansion 



f 2 (n 

V'(A) = V'(iPo) - 4>o^V"{^) + ... = -Vo— 



then we have 



V l (a) = -m l %^^ + ... . (16) 



We also need the second derivative of the potential with respect to the order parameter. Taking the derivative of 
Eq.([l5|) with respect to (jj, we obtain 



where the transverse projection operator is defined by 

Pij [rfi] — Sij — rhiirij 
Evaluating Vy [m] in the bulk, where we can use Eq. rtl4) , we obtain 



Vij[m] = qlrhiiJij - P^ — . (17) 

where q 2 = (n — l)/£ 2 . 

Armed with these results, we next discuss the equations of motion governing the fields m and u. The idea is to 
separate the original order-parameter equation of motion into equations for m and u which insure that we do obtain 
ordering, fluctuations are small and the zeros of rfi reflect the zeros of the order parameter ip. The condition that a 
govern the ordering requires that in the bulk, away from any defect cores, t? 2 — ► tpQ and u 2 — ► 0. We expect, for both 
the scalar and continuous cases, that u controls the transition region between the defect core and the ordered bulk. 

Let us look first at the simpler scaler order parameter case. If the form given by Eq. (|10| ) is inserted into the equation 
of motion given by Eq.(0), we obtain, without approximation, the equation of motion for u: 

A(l)u(l) + <j 1 (l)A(l)m(l) = -V'[a(l) + u(l)}+a 2 {l)(Wm{l)) 2 . (18) 

We assume that the equation of motion satisfied by m is of the form, 

A(l)m(l) = 5(1) , (19) 

where 5(1) is a functional of m(l) and, if naive scaling is to hold, we require S w Q(L~ l ). How one proceeds when 
one does not have naive scaling will be discussed elsewhere. Using Eq.(n3) in Eq.dia) for u leads to an equation of 
motion for the field u(l): 

A(l) U (l) = -y'[a(l)+ U (l)]+a 2 (l)(Vm(l)) 2 -a 1 (l)S(l) . (20) 

Our goal then is to show that we can choose 5(1) such that u is small in the bulk and near a defect. More 
quantitatively, in the bulk, we have a[m] — * ^osgn[m], while lim^^^ulm] = 0. To see how this works and to put 
some constraints on 5(1), we note, if u is small in the bulk, then we can self-consistently expand 

V'[a(l) + u(l)] = V'[a(l)} + V"[a(l)]u(l) +... = a 2 (l) + q 2 u(l) + ... . 

The equation of motion for u then takes the form in the bulk: 

[A(l) + q 2 } «(1) = -,7 2 (1) [1 - (Vm(l)) 2 ] - o x (1)5(1) . (21) 

In working in the bulk ordered regime, which makes the dominate contribution to the scaling properties, we can 
estimates m « L,V « L~ x , and d/dt w L~ 2 . Notice on the left-hand side of Eq.(pi|) that u has acquired a mass and, 
in the long-time long-distance limit, the term where u is multiplied by a constant dominates the derivative terms and 
u is given by 

q 2 u{\) = -a 2 (l) [1 - (Vm(l)) 2 ] - ai(l)S(l) . 

In the limit of large \m\ the derivatives of a go exponentially to zero and the right-hand side of Eq.(E^) is exponentially 
small. Clearly we can construct a solution for u which is also exponentially small in the bulk. We have then on rather 



general principles that the held u must vanish rapidly as one moves into the bulk away from interfaces. We then have 
the rather weak constraints on 5: 

i. 5 must be chosen such that to grows and the held 3 orders. 

ii. If the system satisfies naive scaling then 2 must go as C(L _1 ) in the bulk. 

The form for 2 in the bulk which fulfills these requirements is given by 

5(1) - sgn(m(l)) (g (l) + 5l (l)(Vm(l)) 2 + ■ ■ •) 

where <?o(l) and <7i(l) must go as L~ x for long times. 

Let us now look at how this separation process carries over to the case of a continuous order parameter n > 1. 
Inserting Eq.(nOh into Eq.(pJ) generates an equation of motion for u(l) of the form 



A(l)«i(l) = -Vi(ff(l) + u(l)) - A(l)cr<(l)] = -Vi(a(l) + 5(1)) - a i;j A(l) mj (l) + a % , ]k V a m 3 V a m k . (22) 

We assume that to satisfies the nonlinear equation of motion 

A(l)m(l)=S(l) , (23) 



where 5 has the same interpretation as in the scalar case. Inserting Eq. (E3J) back into Eq. (E2|) gives the basic equation 
of motion for u : 

A(l)Uj(l) = -V r 1 ((7(l)+M(l))-CT l;J -S J (l)+(Tj ;jfc V«TO J -V a TO fe . (24) 

In the bulk we self-consistently assume that the potential can be expanded in powers of u: 

V i (a(l)+ti{l)) = V i (a{l)) + V ij {a(l))u j {l))+... . (25) 

The two terms on the right-hand side of Eq.(|25|) can be simplified using Eqs.([lq) and (|lj) and lead to the new form 
for Eq.(||): 

A(l)«i(l) = g ™» _ <?0™i u i "I 2 U ^ _ ^J^iC 1 ) + CTiyTcVaTO^VaTOfc (26) 

where we have divided u into its longitudinal and transverse parts: 

m = rhiUL + uj , (27) 

where to • u T = 0. Dotting rhi into Eq.(|2q), we obtain the equation determining the longitudinal part of u, 

(n — 1) _ _ 

rhiA(l)ui(l) = — 32 Qq u l + mj0iyfcV a TOjV a mfc - mi<j i; jEj . (28) 



m 
In the bulk, starting with ipi = iporhi, we easily obtain 



^° d r-i 
a i;j = — p iA m \ 
m 



rhi&vj ~ 0(L 3 ), mjCTfyHj ~ 0{L 4 ) and term proportional to rhiOij can be dropped in Eq.(28) compared to the 
terms of C(1/to 2 ). Next, since 

Vi;jk = 2 [Sii^k + Siklflj + SjkThi ~ Z^rhjThk) , (29) 

in the bulk, we have 

TOjCrjjjfcVaTOoVaTOfc = ^Pj k V a TO,- S7 a m k . 

TO 2 

Using this result back in Eq.(p8|), the equation of motion in the bulk for the longitudinal part of the fluctuation field 
is given then by 

miA(l)ui(l) = ^ qlu L ^P jk V a mjV a m k ■ (30) 

m 2 m 2 



All derivative terms in Eq.(30) acting on u can be dropped in comparison with the q 2 u^ term in the scaling regime, 



and, to lowest order in L , we can express ut explicitly in terms of m 

1 



%u L = — [(»- 1) — PjkV a mj V a m t ] . (31) 



Because of the mass term, q$ > 0, we see indeed that ul is of order L 2 and no additional constraints are put on E. 
Let us turn next to the transverse part of u. Multiplying Eq.(£9|) by the transverse projector P we obtain 

(n — 1) 
PijA(l)uj(l) = j — uj - PifeUfcjSj(l) + Pi£ae-jkV a mjX7 a m k ■ (32) 



Then in the bulk 



and, using Eq. 



P ik * kj E j (l) = '^P ij E j (l) = ^(l) , (33) 

to to 



ifjQ 

PtfVtukV a m,jV a m k = -2^(m k \/ a m k )P ij V a m J . (34) 



Putting Eqs.(||) and @ back into Eq.(||) gives 

P y A(l) Uj (l) = ^—f^-uj - — Sf (1) - 2^{m k V a m k )P ij V a m j 

which governs the transverse fluctuations. For self-consistency u T must be small in the bulk. This requires us to 
choose the transverse component of S to be given by 



Sj(l) - 2—(m k V a m k )P lj V a m j = 



in 



2 

Sf (1) = (TO fc V Q TO fc )PyV Q mj . (35) 

m 



With this choice, the equation for the transverse fluctuations is given by 

P 4i A(l) Ui (l) = ^P-uJ . 
Thus we have that uj is generated by any coupling back to ul « 0(1/ L 2 ) via PijA(l)uj(l). We can estimate 



.t 



U; 



L 2 ^-A(l)u L (l) « 0(1/L 2 ) , 



and generally u sw C(l/L 2 ). The requirement that the bulk part of the transverse fluctuations u T be small fixes the 
form of S T to be given Eq.(p5|). This form does not depend on any details of the potential, and can be simplified. 
Consider 

2 12 

—m k V a m k = -^-V q to 2 = — V Q m 
to m z m 

and 

<9to 

h 

Inserting these last two results back into Eq.(B5|) gives 

2 

2 i (1) = (V a m)mV Q mj = — 2V Q mV Q mi . (36) 

to 



PijV a rrij = rn— Vq-TOj = m\7 a rhi 



The equation of motion satisfied by m is given by Eq.(|23|). While the transverse part of 5 is given by Eq.fl3q), the 
longitudinal part of 3 is constrained only by the requirement that it scale as L~ l . The precise form for 3/, in the 
scaling regime must be determined self-consistently within perturbation theory. If we look at the building blocks in 
the problem we see that the quantities which are of 0(1) are rh and V a m,i. Thus one sees that the structure of the 
longitudinal part of the equation of motion in the bulk scaling regime can be assumed to be of the general form: 



Sf(l)=7h i (l) 



fl-o(l) + fifi(l) (Vom^l)) 2 + g\) u (l)V a m i (l)V Q m J -(l)V /3 m fc (l)V / 9m < (l) + 



Clearly in the long-time limit we require the g's be of O (L~ x ) and, as we shall see, that we can self-consistently 
construct the <7 p 's if we assume that g p = 0(p) in the vertex expansion. Our final results at second order in our 
expansion will depend on go and g\ . 

The assumption we make here is that the higher-order terms proportional to </' ', for t > 0, contribute in a non- 
trivial way starting at order t + 2. Thus g^ 2 ' , due to various contractions, acts at second order only to renormalizc 
So- 

IV. FIELD THEORY FOR AUXILIARY FIELD 

The equation of motion satisfied by the ordering field rh including terms which contribute up to second order is 
given in the bulk by 

A(l)mi(l) = go(ti)mi(l)) + 3,(1) (37) 

where 3,(1) = Sf (1) + Sf (1), and 

Sf(l)=.gi(ii)m 4 (l)(V Q m,(l)) 2 

Ej (1) = -2V Q m(l)V Q m l (l) . 

The functions <7o(^i-^| an d 3i(*i) are determined within perturbation theory. Our analysis will foUpw the standard 
Martin-Siggia-RoseE3 method in its functional integral form as developed by DeDominicis and Pclitit-3 and presented 
in detail in I. In the MSR method the field theoretical development requires a doubling of fields to include the response 
field M. As in I, we introduce a field h(l) conjugate to m(l) and a field H (1) conjugate to M(l). 

Following closely the formal development in I, we find that the fundamental equation satisfied by the average of 
the field m, in the presense of sources, is given by 

i [A(l) < m(l) > h -Q(l)] = - I d2 n (12) < M(2) > h +H(1) (38) 

where the vector labels are suppressed, 

n#'(12) = 5(h - t )S(h - t 2 )g(n - r 2 )S tJ . 
It is assumed here that the field m(l), at the initial time to has gaussian statistics with variance 

< m Q (fi)m (r 2 ) >= Sijg(n - r 2 ) • 
The nonlinear vertices in Eq. (p8|) are given by 

Qi(l) =< 3,(1) >^ Qf (1) + Qf (1) + Qf (1) 
with 

Qf (1) =< Sf (1) >= go(l) < mi(l) > 

Qf (1) =< 5 4 L (1) >= , 9l (l) < mil) (Vm,(l)) 2 > 



(39) 



QT(1) =< -f (1) >= (-2) < V a m(l)V„»hi(l) > . 
The fundamental equation satisfied by the average of the MSR response field is given by 

- i [A(l) < M(l) > fc +Q(1)] =fc(l) 

where, we define 

and the nonlinear contributions are given by 

Qi(l) =< 5,(1) >= Qf (1) + Of (1) + Of (1) 
with 

Of (1) =< H? (1) >= fld(l) < M^ljpijfl) > 

Qf (1) =< Sf (1) >=< Mj(l) 5l (l)py(l) (V Q m fc (l)) 2 > - < V Q [Af,-(l) 5 i(l)rh i (l)2V a m,-(l)] > 
and 

Of (1) =< Sf (1) >=< 2m,(l)V Q [M J (l)V Q m J (l)] + 2py(l)V a [Mf(l)V a m(l)] > , 
and 

Ptj'C 1 ) = — (^y ~ Wjmj) . 

All correlation functions of interest can be generated as functional derivatives of < m(l) >;,. or < M(l) >/j with 
respect to h(l) and H(l). In the limit in which the source fields vanish, each term in the two fundamental equations 
vanish. Therefore it is derivatives with respect to the external sources of these equations which are of interest. Let us 

introduce the notation that Gax.a 2 ,A„(12...n) is the n th order cumulant for the set of fields {Ai, A2, ...., A n } where 

field Ai has argument (1), field Ai has argument (2), et cetera. This notation is needed when have cumlants with m 
and M mixed. As an example 



G 



Mi, 



,(1234) 



S 3 < m(4) > h 



SH(l)Sh(2)Sh(3) 
As a short hand for cumulants involving only m fields we write 



G (n) (12- 



8h(n)5h(n- I) ■ ■ ■ 5h(2) 



< m(l) >h 



(40) 



(41) 



The hierarchy of equations connecting these cumulants is given by taking functional derivatives of the fundamental 
equations given by Eqs.(pq) and (39). 

The equations governing the n'^irder cumulants are given by 



A(1)G 



Mi, 



(12...n) + Q n (12...n) 







and 



i A(l)G (n) (12...n)-On(12...n) 
where the Q's are defined by 

0n(12...n) 



dl n (il)G Mm . 



,(12.. .n) 



rn — 1 



5h(n)Sh(n - 1) • • • 6h(2) 



0(1) 



(42) 



(43) 



(44) 



and 

Q "( 12 --^ ^M^-I)...^(2) Q(1) ' (45) 

With this notation the equations determining the two-point functions can be written as 

- * [A(l)G Mm (12) + Q 2 (12)] = 5(12) (46) 

t [A(1)G(12) - Q 2 (12)] = - /rfin (ll)G Mro (12) (47) 

The point now is to show that there is a consistent perturbation expansion for this theory. To get started we need 
to express Qi(l) and Qi(l) in terms of a fundamental set of vertices which can be written in terms of the singlet 
probability distribution 

P h (x,l)=<6(x-m(l))> h . 

After a great deal of rearrangement one can show that the nonlinear vertices, the Q's can be put in the form 

QP(l)=go(l)Ui(l) (48) 



Qf(l) = J2 [ dm 01(1)^(123)0^(23)^(1) 

3 J 

Qf(i) - J2 I dm w(m)o lk {ii)u ]Kl (i) . 

jk 

Qf( 1 ) = -E/^3go(l)0^.(l)C/ 4J (l) • 

3 

Qt'H 1 ) = -J2 / d ^3 ff i(l)«;(123)Ofcfc(23)0^(l)^(l) 

jk 

Qf^i 1 ) = ~J2 / d2 ~ d3 2g 1 (l)w(123)O i (2)0 Mj (S)U j (3) 

3 

QVi 1 ) = -J2 I dm w(123)0 Mj (l)O se (23)Q siiij (l) 

3 si J 

Q^i 1 ) = -J2 I dm 2w(123)Om,(3)0^(2)[/^(1) (49) 



j S e 

where we have introduced the operators 

5 
5h~(2) 

and 



^ = J^+ Gi 3 l) V) > 



°^ (23) = S hj (2)6h k (3) + <%<& + G ? } (2)^5) 



<Vo\ 6 i /n,(i),o^(i), 



We have also introduced the three-point vertices 

d 

W (123) = ^VW<5(12)VW^(13) 

Q = l 

and 

u)(123) = V tt (1) [(5(13)V Q (1) ^(12) 

Each term in these expressions for the Q's can be expressed in terms of the set of nonlinear vertices which are integral 
moments of P^ [x] : 



U ijk... ,£mn 

We have also defined 



...(1) = I d n x x lXj x k . . . V*V™V£ . ..FWS, 1] . (50) 



Grf,«(l) = U sl>ij (l) - U is>ij (l) - U u , sj (l) 



V. PERTURBATION THEORY EXPANSION 

All of the cumulants involving the field rh can, in principle, be obtained from Eqs.pq) and ( p9[ ) by taking functional 
derivatives. This then requires that we work out the functional derivatives of Q n and Q n which are defined by Eqs.(|44|) 
and (Eq). These objects are functional derivatives of Q\ and Q\ which are proportional to a few of the Uijk...,lmn... (1) 
and functional derivatives of these quantities. From this discussion it should be clear that all of the Q n and Q n can 
be written as a product of cumulants multiplying vertices given by 

AAA 

U ijk ... itmn ... iStu ...(l; 234 . . .) = a^(2) ^(3) ^(4) ' ' ' U W-*™-0) ■ 

The point we want to establish is that if U has p vector labels, then, at lowest order, we can take U to be of 0(p/2 — 1), 
plus higher order terms. 

The perturbation theory expansion for the Vijk...,tmn...(X) follows from the expansion properties of the singlet- 
distribution function. The perturbation theory expansion for Ph(x, 1) is straightforward. Using the integral represen- 
tation for the (5-function, we have 

where 7i(l) = ik ■ m(l). The average of the exponential is precisely of the form which can be rewritten in terms of 
cumulants: 



$(£,!) =< e w(1) > h =exp 



OO 1 



s 

8=1 



where G^ (1) is the s th order cumulant for the field H(l). Since Ti(l) is proportional to m(l) these are, up to factors 
of ik to the s th power, just the cumulants for the m field: 
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^ ) (l) = (*) a E*-i^ G 2?a 3 ( 11 ) 



G^(l) = « 3 E ^^^G^^llll) 



and so on. We can therefore write 

$(£, 1) = ezp 



00 (: \ s 

/ , —j-k ai k a2 ■ ■ ■ k a3 G aict2a3 Q , e (ll...l) 



We will assume, as we will show self-consistently, that p order cumulants are of order ^ — 1 in the vertex expansion. 



Expanding and keeping terms up to the 4-point cumulant, we obtain 



Ph(S,l) 



1 E ^£) a2Q3 (m)v^v^v^ 



where 



and 



3! 



p{ n\^)=j^f^y- i% - s 



4 (o) ( j5 )1) = e <E^^(i) e -iE ei ^^GS. J (") 



if(*,l) 



We can define the lowest order set of vertices 



U 



(o) 

ijk...,£mn. 



(\) = J d n X XiXjXk . . . ViV™V£ ■ ■ ■ P { h ] % 1] 



(51) 



(52) 



d n xx l x 3 x k ...ViV™V™. 



(2tt)™ h V ' ' 



\Jb Ju iL> <t JLj 7 X £* 



(2tt)" 



(-iki)(-ik m )(-ik n ) . 



d n k 

(2tt)™ 



i(°) 



'/k-x 



^(*,l)e-* 



It should be clear, after inserting Eq.(pl|) back into Eq.(pfj|), that 



O1O2O3O4 



ijk...,lmn 



a) 1 v^ r/(°) n^— <7^ (11111 + 

...V / ' / ,, ij'fc. .. ,£mn. .. 01O2O3O4 V x 7 ^| v - T oio 2 0304 v^- 1 -- 1 -- 1 -/ ~ 



(53) 



(54) 



It is clear that if the term on the first line of Eq.(p4|) is of 0(p), then the term on the second line is of 0(p + 3) 

For these ideas to be self-consistent then the functional derivatives of U^ . . . must be higher order. The reason this 
works is because factors of the one-point cumulant Gi(l) enters in the exponential appearing in $i (fc, 1) multiplied 

by a factor of fc. Thus functional derivatives either bring down factors of k from the exponential or increase the order 
of cumulants which do not involve the one-point cumulant. To see how this works consider the set of derivatives 



C/... Q3Q4 (1;34) 



■M, 



6h a3 (3)6h a< (4) 



tf (0) (l) 



The functional derivatives then act on $^ (k, 1). It is then easy enough to work out, using Eqs.(|53|) and (p2J), that 
6 2 



Sh a3 (3)Sh ai (4) 



4 0) (£, 1)U =0 = E k ^ k ^ 0^(13)0^(14) - ^^^(1134) 



$$>\k,l) 
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We then have 



U 



(1;34) = J2 U {0 l ia2 (l) ^3(13)0^(14) - ±G£> QaQ3Q4 (1134) 



Q1Q2 



and if C/. (0) (l) is of (9(p) then EA u i ia2 is of (9(p + 1), G^( a2a3CLi is of 0(1), and U.. a3(Xi (l; 34) is of 0(p + 1) plus 
higher order terms. 

Notice then that all derivatives of U and higher order terms can be expressed as products of the vertices C/' ' 
with legs connected by cumulants. Any contribution to the nonlinear vertices Q n and Q n can therefore be ordered 
summing up the contributions from each vertex and cumulant where a vertex with 2s + 2 legs makes a contribution 
of O(s) and a cumulant G ( - 2p+2 - ) makes a contribution of 0(p). 

For our purposes here, we only need two sets of the bare vertices U^ ': 



V s (ai;a 2 ,. . .a 2s+2 ;l) = / d n x x ai 



d n k 



K i)k a2 k a3 . . . k a2B+2 'P h {k,l)e 



-ik-x 



and 



V s (ai,a 2 ; ...a 2s +2;l) = / d n x x ai x a2 / k a . 3 k a4 . . . k a2s+2 <&^> (k, l)e lfc ' x . 

It can be seen the V s (ai;a 2 . . .a 2s + 2 ; 1) is the natural generalization of the vertex </> s (l) given by Eq.(92) in I. As 
in I, it is straightforward to work out these vertices explicitly by doing the integrations. The quantities we will need 
below are given by: 



Vo(ai;a2; 1) = $ ai ,a 2 



I 2 r(gffi 

5(1) nr(f ) ' 



Vi(ai;a 2 ,a 3 ,a4; 1) = I ai ,a 2 ,c 



I 2 r(ffi) 

S 3 (l)n(rH-2)r(|) ' 



and 



V (ai,a 2 ; 1) = S aia2 



Vi(ai,a 2 ;a3,a 4 ; 1) 



S(l) 



Oai,a 2 Oai,a 2 , rt -Lai,a 2 ,a 3 ,aA 

n ~ n + 2 ' 



W have also introduced 

S(l)^G^(ll)=<mi(l)> 
with no summation over a, and the symmetric tensor 



L ai ,a 2 ,a 3 ,ct4 



0ai 1 a 2 0a 3 ,a4 t ®ol\. a 3 ^a 2 . ct& T" ai ^ a4 a 



A. Two-Point nonlinear vertices to second order 



Armed with these results we can return to the evaluation of the Q n and Q n derived from Q\ and Q\ which are 
given by Eqs.([48|)-(p}). We summarize first the results for the two-point quantities which enter Eqs.(E6|) and(E^). We 
find after a significant amount of work that 



Q 2 (12) = ft(ll)G(12) + S*(1234)G(2342) 



(55) 



Q 2 (12) = fi(ll)G M m(12) + 5 1 (1234)G Mmmm (2342) 



(56) 
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where integration and summation over repeated barred indices is implied. The various auxiliary quantities are defined 
by: 

0(12) = ( 50 (l) + <?i(l)S (2) (l)) V (ai;a 2 ; 1)5(12) (57) 

5(1234) = ^V (1234) + Vl(1234) + Vi (1234) (58) 

Si (1234) = ^Vb(1234) + Vl(1234) + V 2 (1234) + Vs(1432) + V£ (1432) (59) 

Vb(1234) = -$t,(l)Vi(ai; a 2) a 3 , a 4 ; 1)5(12)<5(13)<5(14) 

14(1234) = .gi(l)Vo(ai;a 2 ;l)5 Q3Q4 ^(134)5(12) 

Vi(1234) = Vi(a3, 04; oi, as; l)5 a3ai W {134)8(12) 

V 2 (1234) = -W(lM)5(12)QW teUtaitaa (l) 

1/ 3 (1234) = -2Vi(ai,a 3 ;a 4 ,a 2 ; 1)1^(134)5(12) 
Vf (1234) = -2 5 i(4)Vo(a 4 ;a 2 ;4)<J Ql , Q 3 1^(134)5(24) 



and 



2L 3 a 4 ,a 1 ,Q :2 ( 1 ) = -Vi(a 3 , «4! ai, a 2 ; 1) + Vi(ai, a 3 ; a 4 , a 2 ; 1) + Vi(ai, a 4 ; a 3 , a 2 ; 1) 

<5 (2) (1) = ^ < V l m Q (l)V,m a (l) > 

i 

where there is no summation over a. These results give a closed solution for the two-point correlation functions at 
zeroth-order. This analysis of this lowest-order solution is given in the next section. In the following section we 
analyze the four-point correlation functions needed in order to extract the two-point correlation functions at second 
order. 

VI. ZEROTH ORDER THEORY FOR TWO POINT CORRELATION FUNCTIONS 

The equations of motion at zeroth order for the two-point correlation functions are given by the coupled set of 
Eqs.fliq) and ( (47|) with the lowest order contributions for Q 2 and Q 2 given the leading order terms in Eqs.(J55|) and 
(pq). Inserting 

Q^ 0) (12) = f7(lI)G(12) 

Q^ 0) (12) = O(lI)G M m(12) 
into Eqs.(f46|) and ( j47| ) and explicitly writing the vector labels, we obtain 

- i [A(l) + wo(l)] G { °l ima2 (12) = 5(12)5 Ql , Q2 
13 



i [A(l) - «o(l)] Gi°) Q2 (12) = -J dl n (ll)G^ ima2 (12) , 

where we have defined 

^o(l)=(5o(l)+3i(l)5 (2) (l))v (a 1 ,« 1 ;l) . 
We see at once that the solution to this set of equations is diagonal in the vector-indices, 

and 

^ r M c « 1 m„ 2 (12) = oli, oc 2 G Mm {12) , 

where G < -°- ) (12) and G^ fm (12) are the same quantities found in the scalar case in I. We summarize briefly the results 
since they are needed here. At long times we can write 

Mt) - ^ (60) 

where u is a constant we will determine and t c is a short-time cut-off which depends on details of the early-time 
evolution. One has then that the response function is given by 



G^M^HG^MAH-^-O'^^^ ' 



where r — f± — r*2, the correlation function is given by 

where g(0) is the on-site value of the initial correlation function, and it is convenient to define 

T= h+t 2 
2 

If we are to have a self-consistent scaling equation then the autocorrelation function (r = 0), at large equal times 
t\ = t<i = t, must satisfy 

2u 



**<«> -««•>(£) (Sf s * ■ 



Clearly this result fixes the exponent 



and the amplitude 



I, d N 

^2 (1+ 2 } 



1 ff(0) , fi9 . 

= (i c ) 2 - WF 2 ' ( } 



Eq.(pl|) can be rewritten in the convenient form 

Gt 0) (r,ti*2) = V /5W(ti)5(o)(t 2 )$ ( o ) (ii,t 2 )e-' r2/(4T) 
where 

d/2 



'I'imiI'i-':'! = ! ^7^ 
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The nonequilibrium exponent is defined in the long-time limit by 

G<°>(0,ii,t a ) * ,, , v /V*Ife 



*(o)(*ij*a) 



A 



VS"(«i)50(t 2 ) w V T 

and we obtain the OJK result A = |. Looking at equal times we have the auxiliary field scaling function 

/o(x) = G( ° )(r ' tt) = e-^ 

where the scaled length is defined by x = f/L(t), and the growth law is given by L 2 (t) = At. The exponent u, defined 
by Eq.fln]), is zero in the OJK approximation. 

A. Four- Point Correlation Functions at First Order 

If we are to evaluate G^ and Gum to second order, we see that we must evaluate the four-point quantities G^ 
and GMmmm to first order in the vertex expansion. This requires the evaluation of Qi and Q4. Using the same 
techniques developed in evaluating Qi and Qi we find 

Q 4 (1234) = fi(ll)G M mmm(1234) + 5i(1234)P M (234,234) (63) 



where £7(12) and 5*1(1234) are defined by Eqs.(57) and (K9h respectively, while 



P M (2'3'4', 234) = G Mm (4'2) [G(3'3)G(2'4) + G(3'4)G(2'3)] 

+ G Mm (4'3) [G(3'2)G(2'4) + G(3'4)G(2'2)] 

+ G Mm (4'2) [G(3'2)G(2'3) + G(3'3)G(2'2)] . (64) 

We also have 

Q 4 (1234) = f7(ll)G (4) (1234) + 5(1234)P(234, 234) (65) 

where 5(1234) is given by Eq.(pq) and 

P(2'3'4', 234) = G(4'2) [G(3'3)G(2'4) + G(3'4)G(2'3)] 

+ G(4'3) [G(3'2)G(2'4) + G(3'4)G(2'2)] 

+ G(4'2) [G(3'2)G(2'3) + G(3'3)G(2'2)] . 

Inserting Eq.(p3|) for Q4 into Eq.([42|) we see that we can do a partial integration and write 

G Mmm m(1234) = G { ^ m {ll)iS x (1234)P M (234, 234) (66) 

where S\ is defined by Eq.(pq). Using the symmetry properties of Pm (234, 234) we can show that Eq.(|6q) can be 
written as 



where 



G Mm mm(1234) = Gf Im {ll)V s (2- 134)P M (234, 234) (67) 



V s {2- 134) = ^Uo(1234) + iV L , s % 134) + iV hs (2, 134) 



and the symmetrized vertices are given by 
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V L<8 (2, 134) = y L (2134) + 14(2, 413) + Vl(2,314) 

and 

Vi )S (2,134) =Fi(2134) + Fi(2,413) + Fi(2,314) . 

In turn Eqs.(|67|) and ( |65| ) can be put back into Eq.p^) with n = 4 to obtain G^ . After manipulations it can be 
written in the properly symmetric form 

G< 4 )(1234) = iG^(ll)K(l;234)P(234,234) + ±0^(21)^.(1; 234)P(234, 134) 

+ ^ } M (31)K(1; 234)P(234, 124) + ±G^(41)V S (1; 234)P(234, 123) . 
We see after these manipulations that all of the first order corrections have been combined into a single vertex. 

VII. TWO-POINT CORRELATION FUNCTION AT SECOND ORDER 

A. General Equations 

Given G^ and GMmmm at first order, we can return to Eqs.(|63|) and ( |56| ) to obtain Q 2 and Q2 to second order. 
These in turn are put back into Eqs.(p6[) and (|47]) to obtain the second-order results for the two-point correlation 
functions. We focus here on the correlation function. After a single integration of Eq.([47|) we have 

G(12) = G^(12) + G%(li)^V s (l; 234)G 4 (2342) + G^(11)V S (2; 134)G M mmm(2342) 
= G(°)(12)+G^(11)^K(1;234)^ S (1';2'3'4') 

G^ M (21')P(2'3'4', 234) + G% (2l')P(2'3T, 234) + G^ M (3T)P(2'3'4', 224) + G^(4l')P(2'3'4', 223) 

+ G(°)(11)K(2; 13l)G^ m (2l')V s (2'; l'3'4')P M (2'3'4', 234) . 

This last term simplifies since, because of causality, only the term proportional to GMm{2'2) in Pm(2'3'4', 234) 
survives, and 

V S (V; 2'3'4')P(2'3'4', 234) = 6V S {V; 2'3'4')G(2'2)G(3'3)G(4'4) . 

We then have the final formal expressions: 

G(12) = G (0) (12) + G (s) (12) + G {u \l2) + G [u) (2l) 

where the symmetric contributions are given by 

G {s \12) = \g%{11)G ( ^ m {21')V s {1-WA)V s {1'-2'^ 

and the unsymmetric contributions are given by 

G( u )(12) = 2G|"] Vf (ll)G( H220T4(l;234)T/ s (1^2 , 3'40G^] v/ (2lOG( )(3 , 3)G^(4'4) . 

The detailed analysis of these contributions to the correlation function follows closely the analysis developed in 
detail in I. Indeed if the full vertex V 8 is replaced by (i/2)Vo and n set to 1, these equations reduce to those found 
in I. One can again carry out explicitly the internal spatial integrations. Among the new elements in this analysis is 
the treatment of the gradient insertions in the vertices, and the internal vector sums. One must also introduce the 
parameter, 

ff = 5i(l)Vo(ai;ai: 1)5(1) 

which, as anticipated earlier, requires that 31 (1) go as L~ l for long times with an amplitude which is determined as 
part of the scaling structure. 

1G 



B. Extraction of Indices 

As in I, all of the various log singularities found in second order, and arising from internal time integrations, can 
be absorbed into expressions for the exponents u>, A and v. 

At second order in the vertex expansion the exponents are determined by the set of equations: 



x d 2 2 d+1 tJ 



and 



= co 2 2 d+1 



Hu 



3 d / 2 



(68) 



(69) 



where the quantities Hjj and Hs are given below. The condition that the growth law L sa t 1 ' 2 be maintained order 
by order in perturbation theory determines the parameter lo, and, as in I, can be expressed in terms of the exponent 



2u) + - = 1 + d/2 



(70) 



There are two intrinsically different contributions, Hjj and Hs, to the second-order expressions for the correlation 
functions which come from graphs with different structures. Hjj and Hs, which depend only on the parameters u>, g, 
d and n, are defined in terms of a set of auxiliary quantities: 

H s = QfM d + QfM d+2 



Hu = Q^K d + QfK^ + QfK^ 



where 



Qf 



2(n 



+ 2)\ oj J 



)(2) 



18lu 2 



(d + 5)S (3) +2(d-l)S {4 '> 



qP = 



2(n + 2) V w 



1 - *S™ 



?£> = 4^ (3) 



(9) - — 



S& 



and the S^ % ' are basically the result of internal n-vector sums: 

S (2)_ 2(n-l) | (n + 2) ( 
in 3 



S (3) = 2(n-l) 
n(n + 2) 



ng 



g (4)_ (n-2)(n-l) 2(n-l) 2 

n^(n + 2) n 
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The constants g and ui, parameterize the scaling properties of the nonlinear terms in the equation of motion for the 
auxiliary field. Finally we have the <i-dependent integrals 



Ma 



dz 



z 4 ' 2 - 1 ir 2 (rf/2) 

[l + z] d ~ 2 T(d) 



(71) 



K 



(o) 



dz 



y d/2-l 



[(I + z)(3 - z)} d / 2 



(72) 



K 



,(6) 



dz 



y d/2-l 



[(i + z)(3 - z)] d / 2 (i + z y 



2(1 



(3-z) 



(d + 2 



{1-zf 



{i-zf 



K (9) = _d 



dz 



z 



d/2-l 



4j [(1 + *)(3 -*)]*/» (l + z)(3-*) 



l-(d + 2) 



(1-*) 

(3-*) 



(73) 



If we set n = 1 and 5 = then this set of equations reduces to that found in I. Then Eqs.(p9[) and (f70j) can be solved 
for ui and v and the results inserted in Eq.(p8|) to obtain the index A as given in I. 

The parameters ui and g should be thought of as being generated by some type of renormalization group analysis. 
Carrying forward the RG analogy, these parameters are to be determined as part of finding the scaling fixed point 
in the problem. This process is similar to finding a fixed point Hamiltonian in critical phenomena. The parameter 
u] occurs naturally at lowest order in the perturbation theory expansion, while g naturally arises at second order in 
the expansion, u) and g are determined by the requirements that the scaling law L ss t 1 ' 2 and the index relation 
v = 2A — d hold at all orders. The maintenance of the growth law leads to the condition given by Eq.(|7C|). The 
requirement v = 2A — d is enforced by choosing 



Hu = 



(74) 



While there are many possible ways of extracting explicit numbers for the indices from the perturbation expansion 
just described, we discuss two here. In the expansion method we set u> = ujq — i(l + |) in the second-order terms and 
obtain the indices directly. In the second method we look for a self-consistent solution of Eqs.(pOf) and (I74J) for g and 

u>. For large d and n these two approaches are equivalent. While the various (i-dependent integrals, K d , etc. can be 
worked out analytically for specific values of d, the expressions are not very illuminating. Numerical values for A, v, 
and g are given in Tables /, 77, and III. Except for the values of A maked by an *, whose significance is discussed 
below, the values of v are given by v = 2 A — d and ui is given by Eq. (|74|) . One sees that the self-consistent values for 
A are all close to the OJK values. The perturbation theory results can lead to much larger corrections. 

It is instructive to work out the large-<i limit analytically. For general n, one finds a solution of Eq.(f3) in the limit 
with 
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(n + 2) 



1 2(n-l) 

- + — 

4 3n 



This expression for g has a minimum value of 1/4 for n = 1, a maximum of 7/16 for n = 2 and then a slow decay to 
zero as a function of n. For the scalar case the contribution to v and A = d/2 + v/2, are given to leading order by 

V2^d 3 / 2 



96 3 d / 2 



which gives exponential decay to zero for large d. 

For large n the exponents are also given by the OJK result, with corrections of the form A = d/2 + \i(d)/n + • • •, 
where the precise dependence of Ai on d is complicated. 

This procedure for fixing the coefficients u> and g works straightforwardly for the scalar case and generally for d > n. 
However, for d < n one finds, for small enough d 7 that the ne w sp in-wave contributions (proportional to (n — 1) in 
the £M) lead to a breakdown in this process. Solutions to Eq.(|74|) do not exist and one can not enforce the relation 
v = 2A — d. In this case we have chosen g such that Hjj is a minimum. The structure of the theory for n > d needs 
further work. This is just the regime where one does not generate stable topological defects. 
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dimension 




n= 1 


n = 2 


n — 3 




per 


0.5819 


0.8120 


1.1172* 




sel 


0.5154 


0.5590 


0.6206* 


1 


TUG 


1.0 


0.699 


0.622 




OJK 


0.5 


0.5 


0.5 




num 






0.648 6 




per 


1.0530 


1.2045 


1.5326* 




sel 


1.0059 


1.0227 


1.0597* 


2 


TUG 


1.2887 


1.171 


1.117 




OJK 


1.0 


1.0 


1.0 




num 


1.246 ± 0.02 a 








per 


1.5375 


1.6284 


1.8240 




sel 


1.5024 


1.5082 


1.5216 


3 


TUG 


1.6726 


1.618 


1.587 




OJK 


1.5 


1.5 


1.5 




num 


1.838 ±0.2 a 






large 




d/2 


d/2 


d/2 



TABLE I. Values of exponent A. In second column per refers to values from the current theory fully expanded, sel refers to 
a self-consistent solutions from the current theory, TUG refers to values from Ref.( Ij), OJK refers to the values from Ref.( H), 
and num refers numerically determined values: a is from Ref.( |l0| ), b is from Ref.( H3). An asterisk indicates that no solution 
to Eq.(^4J) was found 



dimension 


n = 1 




n = 2 


n = 3 


1 

2 

3 

large 


7.6423632242e - 
5.1175087572e - 
4.2674144077e - 
0.25 


-01 
-01 
-01 


1.2448881672e+00 

9.3442756917e - 01 

7.9300891113e -01 

0.4375 


1.3940000000e + 00* 

1.2310000000e + 00* 

1.0569791347e + 00 

0.416666 



TABLE II. Values of parameter g 



dimension 


71 = 1 


n = 2 


n = 3 


1 
2 
3 


3.08732e - 02 
1.17040e-02 
4.80000e - 03 


1.181836e-01 
4.530440e - 02 
1.645200e - 02 


1.7438e + 00* 
1.3928e + 00* 
4.3149e - 02 



TABLE III. Values of exponent v 



19 



VIII. CONCLUSIONS 

It has been shown how one can extend the method developed previously for a scalar order parameter to the case of 
the n- vector model. The approach developed here appears to be a rather general tool for looking at field theories where 
the field is growing and showing scaling behavior. One is able to develop a systematic expansion in the number of labels 
on the non-linear vertices appearing in the problem. This expansion leads directly to expressions for the anomalous 
dimensions in the problem. Less generally one is then confronted with the interpretation of the perturbation theory 
expansion in a particular realization of the theory. As organized here, the self-consistent corrections to the OJK 
results for the indices A and v are typically quite small and vanish for both large d and n. The large d convergence is 
tied to the enforcement of the equation v — 2A — d relating the indices. 

The transverse degrees of freedom enter quite differently into the problem compared with the longitudinal degrees 
of freedom. The longitudinal contributions to the non-linear terms in the equation of motion for the auxiliary field 
m must be determined self-consistently in constructing the scaling properties in the problem. The contributions to 
the transverse part of the equation of motion for the auxiliary field are, because the transverse degrees of freedom 
are massless, fixed by the requirement, for self-consistency, that the amplitude for the transverse order parameter 
fluctuations be small compared to the ordered component. It turns out that the transverse contributions to the 
equation of motion for the auxiliary field are sufficiently strong, for fixed n > 3 and sufficiently small d, that we are 
unable to enforce the condition v = 2A — d. This regime requires further study. 

The point of view developed here is somewhat unanticipated. In the most direct approach, as discussed in some 
detail in I, one makes the substitution a = a[rn] into the order parameter equation of motion to obtain the equation of 
motion for m. The path taken here is quite different since the equation of motion for m is constructed self-consistently. 
The surprising point is that the quantity Si(m), which enters the equation of motion satisfied by the auxiliary field, 
is not determined when we insert ip = a + u into the order parameter equation of motion. It is this freedom which 
allows us to construct the scaling regime form for 5^ (fa) . 
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